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ABSTRACT 

We present a homogeneous catalogue for blue straggler, red giant branch, horizontal 
branch and main-sequence turn-off stars in a sample of 35 clusters taken from the ACS 
Survey for Globular Clusters. As a result of the superior photometry and relatively 
large field of view offered by the ACS data, th is new catalogue is a significant improve- 
ment upon the one presented in lLeigh. Sills fc Kniggc (2007). Using our catalogue, we 
study and compare the radial distributions of the different stellar populations. 

We have confirmed our previous result (jKnigge. Leigh fc Silldii2009.) that there is 
a clear, but sub-linear, correlation between the number of blue stragglers found in the 
cluster core and the total stellar mass contained within it. By considering a larger 
spatial extent than just the core, our results suggest that mass segregation is not the 
dominant effect contributing to the observed sub-linearity. We also investigate the 
radial distributions of the different stellar populations in our sample of clusters. Our 
results are consistent with a linear relationship between the number of stars in these 
populations and the total mass enclosed within the same radius. Therefore, we con- 
clude that the cluster dynamics does not significantly affect the relative distributions 
of these populations in our sample. 
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1 INTRODUCTION 

Colour-magnitude diagrams (CMDs) are one of the most 
important tools available to astronomers for studying stel- 
lar evolution, stellar populations and star clusters. And 
yet, there remain several features found in CMDs whose 
origins are still a mystery. Examples include horizontal 
branch (HB) morphology, the presence of extended hori- 
zontal branch ( EHB) stars, and blue stragglers (BSs) (e.g . 
Sandaed [ToHl : IZinn fc Barnej 1 19961 : iPeterson etafl l2003l : 



Dotter et al]|2010l ). Previous studies have shown that the 



observed differences in the HBs of Milky Way globular clus- 
ters (GCs) are related to metallicity ( Sandage & Wallcrstcin 
1 19601 ). however at least one additional parameter is re- 
quired to explain the spread in their colours. Many clus- 
ter properties have been suggested as possible Second and 
Third Parameters, including age, central density and clus- 
ter luminosity, although no definitive candidates have been 
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identified (e.g. iRoodI Il973l : iFusi Pecci et al.lll993l ). An ex- 
planation to account for the existence of BSs has proved 
equally elusive. Many BS formation mechan isms have been 
proposed, including stellar collisions (e.g. [Leonard 19891: 
ISills & Bailvn' '1999^) and binary mass-transfer (|McCreal 
1964; Mathicu & GcUcr 2009). However, no clear evidence 
has yet emerged in favour of a dominant formation mecha- 
nism. 



In short, we still do not understand how many of the 
physical processes operating within star clusters should af- 
fect the appearance of CMDs (e.g. I^ s i Pecci ct al. 19921: 
Buonanno et al.l 1 19971 : iFerraro et al. r Tl999l : |B eccari et al.l 
I2006i ). In general, the importance of these processes can be 
constrained by looking for correlations between particular 
features in CMDs and cluster properties that serve as prox- 
ies for different effects. For example, the central density can 
be used as a rough proxy for the frequency with which close 
dynamical encounters occur. Similarly, the cluster mass can 
be used as a proxy for the rate of two-body relaxation. Once 
the relevant effects are accounted for, CMDs can continue to 
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provide an ideal tool to further our understanding of stellar 
evolution, stellar populations and star clusters. 

It is now clear that an important interplay occurs in 
clusters between stellar dynamics and stellar evolution. For 
example, dynamical models have shown that star clusters 
expand in response to mass-loss driven by stellar evolution, 
particularly during their early evol utionary phases when 
mass i ve stars are still present ( e.g. 'Chcrnoff fc Weinberg) 
ll990l : |Portegies" Zwart et al. I l 19981 : [Heles et al. 2010 '). Mass- 
loss resulting from stellar evolution has also been pro- 
posed to cause horizontal branch stars to exhibit more 
extended radial distributions relative to red giant branch 
and main-sequence turn-off (MSTO) stars in globular 
clusters having short central relaxation times relative to 
the average HB lifetime (e .g. ISigurdsson fc PhinnevI [l995l : 
iLeigh. Sills fc Knig ge'2009'). This can be understood as fol- 
lows. Red giant branch (RGB) stars should be more mass 
segregated than other stellar populations since they are 
among the most massive stars in GCs. HB stars, on the 
other hand, are much less massive since RGB stars experi- 
ence significant mass loss upon evolving into HB stars. Con- 
sequently, two-body relaxation should act to re-distribute 
HB stars to wider orbits wi thin the cluster potent ial relative 
to RGB and MSTO stars (ISpitzer fc Shulllll975l ). provided 
the average HB lifetime is shorter than the central relaxation 
time. Studies have shown that the radial distributions of the 
HB populations in some GCs coul d differ from those of other 
stellar populations. For instance, ISaviane et all l| 19981 ) pre- 
sented evidence that blue HB stars could be more centrally 
concentrated than red HB and sub-giant branc h stars in the 
GC NGC 1851. Conversely. [Cohen et al.1 (|l997l ) showed that 
blue HB stars could be centrally depleted relative to other 
stellar types in the GC NGC 6205. To date, no clear evi- 
dence has been found linking the spatial distributions of HB 
stars to any global cluster properties. 

Peculiar trends have also been reported for the radial 
distributions of RGB stars. For example, a deficiency of 
brig ht red giants has bee n observed in the GC NGC 1851 
(e.g. Ilannicola et al.ll2009l ). ISand guist fc Martell (|2007| ) dis- 
cussed the possibility that this deficiency could be the result 
of strong mass loss on the RGB. Alternatively, some authors 
have suggested that dynamical effects could deplete red gi- 
ants. For instance, giants could experience collisions more 
frequently than other stellar populations due to their larger 
cross-sections for collision (Beers fc Davics 200 4). 

One important example of the interplay that occurs in 
clusters between stellar evolution and stellar dynamics can 
be found in the study of blue stragglers. Found commonly 
in both open and globular clusters, BSs are thought to 
be produced by the addition of hydrogen to the cores of 
low-mass main-sequence (MS) stars, and therefor e appear 
as an extension of the MSTO in cluster CMDs (ISandagd 
Il953l ). This can occur via multiple channels, most of which 
involve the mergers of low-mass MS stars since a significant 
amount of mass is typically required to reproduce the 
obser ved locations of BSs in CMDs (e.g. ISills fc BailvrJ 
Il999t ). Stars in close binaries can merge if enough or- 
bital angular momentum is lost, which can be mediated 
by dynamical interactions with other stars, magnetized 
stellar winds, t idal dissipation or even an outer triple 
companion (e.g. [Leonard fc Linnelllll992l: iLi fc Zhang||2006l; 
iPerets fc Fabrvck^ l2009l : IDervisoglu. Tout fc Ibanoglul 



l201Ci ). Alternatively, MS stars can collide directly, although 
this is also thought to usua ll y be mediated by multiple 
star systems (e.g . iLeonar 3 Il989l : iFregeau et all |2004| : 
iLeigh fc Sillj l20l"ol ). Finally, BSs have also been hypothe- 
sized to form by mass-transfer from an evolving primary 
onto a n ormal MS companion via Roche lobe overflow 
l|McCrea| [r964). 

Whatever the dominant BS formation mechanism(s) 
operating in dense star clusters, it is now thought to some- 
how involve multiple star systems. Th is was shown to be the 
case in even the dense cores of GCs (|Knigge. Leigh fc Silld 
where collis ions between s ingle stars are thought to 
occur frequently ([Leonard ITgsgj ). In [Knigge. Leigh fc Silld 
([20091 ) ■ we showed that the numbers of BSs in the cores of 
a large sample of GCs correlate with the core masses. We 
argued that our results are consistent wit h what is expected 
if BSs are descended from binary stars. [Mathieu fc Gelleil 
([20091 ) also showed that at least 76% of the BSs in the old 
open cluster NGC 188 have binary companions. Although 
the nature of these companions remains unknown, it is clear 
that binaries played a role in the formation of these BSs. 
Dynamical interactions occur frequently enough in dense 
clusters that they are also expected to be at least partly 
responsible for the observed properties of BSs. It follows 
that the current properties of BS populations should reflect 
the dynamical histories of their host clusters. As a result, 
BSs could provide an indirect means of probing the physical 
processes that drive star cluster evolution. 

In this paper, we present a homogeneous catalogue for 
red giant branch, main-sequence turn-off, horizontal branch 
and blue straggler stars in a sample of 35 Milky Way (MW) 
GCs taken from the Advance d Camera for Surveys (ACS) 
Survey for Globular Clusters ([Saraiedini et al.ll2007l ). With 
this catalogue, we investigate two important issues related to 
stellar populations in GCs. First, we test the ob servational 
correl ation found for BSs presented in lKnigge. Leigh fc Sill^ 
(|2009l ) by re-doing the study with newer and more accu- 
rate photometry. The larger spatial coverage considered in 
our new sample offers an important additional constraint 
for the origin of this correlation. Second, we perform the 
same statistical comparison for RGB, HB and MSTO stars 
in order to study their radial distributions. This will allow 
us to test some of the results and hypotheses introduced in 
[Leigh. SiUs fc Kniggel ([20091 ). where we first presented this 
technique for studying stellar populations. In particular, we 
found evidence for a surplus of RGB stars in low-mass GC 
cores relative to MSTO stars. However, we concluded that 
the study needed to be re-done with better photometry. The 
ACS data are of sufficiently high quality to address this is- 
sue. 

In Section (2] we present our selection criteria to deter- 
mine the numbers of BS, RGB, HB and MSTO stars located 
in the central cluster regions. The spatial coverage of the 
photometry extends out to several core radii from the cluster 
centre for most of the clusters in our sample. For these clus- 
ters, we have obtained number counts within several circles 
centre d on the cluster centres provided in Gol dsburv et al.l 
([2010h for various multiples of the core radius. This cat- 
alogue is presented in Section [31 In this section, we also 
present a comparison between the sizes of the different stel- 
lar populations and the total stellar masses contained within 
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each circle and annulus. Finally, we discuss our results for 
both BSs and the other stellar populations in Section [l] 



2 METHOD 

In this section, we present our sample of CMDs and define 
our selection criteria for each of the different stellar popu- 
lations. We also discuss the spatial coverage offered by the 
ACS sample, and describe how we obtain estimates for sev- 
eral different fractions of the total cluster mass from King 
models. 

2.1 The Data 

The data used in this study consists of a sample of 35 
MW GCs taken from t he ACS Survey for Globular Clusters 
jSaraiedini et al.ll2007l )FI The ACS Survey provides unpece- 
dented deep photometry in the F606W (~ V) and F814W 
(~ I) filters that is nearly complete down to 0.2 Mq. In 
other words, the CMDs extend reliably from the HB all the 
way down to about 7 magnitudes below the MSTO. We have 
confirmed that the photometry is nearly complete above at 
least 0.5 magnitudes below the MSTO for every cluster in 
our sample. This was done using t he res ults of artificial star 
tests taken from [Anderson et al.l l|2008l ) , and confirms that 
the photometric quality of the stellar population catalogue 
presented in this paper is very highlj We have also con- 
sidered foreground contamination by field stars, and it is 
negligible. 

Each cluster was centred in the ACS field, which ex- 
tends out to several core radii from the cluster centre in 
most clusters and, in a few cases, beyond even 15 core 
radii. Coordinat e s for the cluster centres were taken from 
iGoldsburv erahl (|2010l ). These authors found their centres 
by fitting a series of ellipses to the density distributions 
within the inner 2' of the cluster centre, and computing a n 
average value. The core radii were taken from lHarrisI l|l99d) . 

2.2 Stellar Population Selection Criteria 

In order to select the number of stars belonging to each 
stellar population, we define a series of lines in the (F606W- 
F814W)-F814W plane that act as boundaries enclosing each 
of the different stellar popul ations. To do this, we fit theoret- 
ical isochrones taken from iDo tt er et all l|2007h to the CMDs 
of every cluster in our sample. Each isochrone was gener- 
ated using the metallicity and age of the cluster, and fit to 
its CMD using the c orresponding d istanc e modulus and ex- 
tinction provided in 'Dotter et alj (|2010l ). The MSTO was 
then defined using our isochrone fits by selecting the bluest 
point along the MS. This acts as our primary point of refer- 
ence for defining the boundaries in the CMD for the different 
stellar populations. Consequently, the selection criteria pro- 
vided in this paper are a sig nificant improvement upon the 
criteria presented in iLeigh. Sills fc Kniggel (|2007l ). and our 

^ The data can be found at 

Ihttp: / /www. astro. ufl.edu/~ata/public_hstgc/ 

Artificial star tests were obtained directly from Ata Sarajedini 
via private communication. 



new catalogue for the different stellar populations is highly 
homogeneous. 

Two additional points of reference must also be defined 
in order for our selection criteria to be applied consistently 
from cluster-to-cluster. First, the selection criteria for the 
HB are determined by fitting a line through the approxi- 
mate mid-point of the points that populate it in the CMD. 
This line is then used to define upper and lower boundaries 
for the HB. Theoretical isochrones become highly uncertain 
at the HB, so it is necessary to specify this additional cri- 
terion by eye. Second, the lower boundary of the RGB is 
defined for each cluster as the point along its isochrone cor- 
responding to a helium core mass of 0.08 Mq. We do not 
include RGB stars brighter than the HB since the tilt of the 
upper RGB varies significantly from cluster-to-cluster, pre- 
senting a considerable challenge for the consistency of our 
selection criteria. Moreover, the distinction between RGB 
and asymptotic giant branch stars in the CMD is often am- 
biguous. 

Example selection criteria for each of the different stel- 
lar populations are shown in Figure[T] Formal definitions for 
the boundaries in the cluster CMD that define the BS, RGB, 
HB and MSTO populations are provided in AppendixlX] We 
note that the sizes of our selection boxes have been chosen 
to accomodate the photometric errors, which contribute to 
broadening the various evolutionary sequences in the CMD. 
With superior photometry, the sizes of our selection boxes 
could therefore be reduced. This would further decrease con- 
tamination from field stars in our samples. 

2.3 Spatial Coverage 

The ACS field of view extends out to several core radii from 
the cluster centre for nearly every cluster in our sample. 
Consequently, we have obtained estimates for the number 
of stars contained within four different c ircles centred on 
the ce ntral cluster coordinates provided in Gold sburv et al] 
We list these numbers only for clusters for which the 
indicated circle is completely sampled by the field of view. 
The radii of the circles were taken to be integer multiples of 
the core radius, and we focus our attention on the inner four 
core radii since the field of view extends beyond this for only 
a handful of the clusters in our sample. An example of this 
is shown in Figure (2] The numbers we list are cumulative, 
so that entries for each circle include stars contained within 
all smaller circles. 



2.4 King Models 

In order to obtain accurate estimates for the total stel- 
lar mass contained within each circle, we generated 
single-mass King mod els calculated using the method of 
ISigurdsson fc PhinnevI (| 19951 ) to obtain luminosity density 
profiles for the majority of the clusters in our sample. 
The pro files were obtained using the concent ration param- 
eters of iMcLaughlin fc van den Marel' j2005l ) and the cen- 
tral luminosity densities of Ha rris (1 9jj3) for each cluster in 
IMcLaughlin fc van den Mare] (|2005l ) that overlaps with our 
sample. We then integrated the derived luminosity density 
profiles numerically in order to estimate the total stellar 
light contained within each circle. After removing clusters 
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Figure 1. Colour-magnitude diagram for the Milky Way globu- 
lar cluster NGC 6205. Boundaries enclosing the parameter space 
in the (F606W-F814W)-F814W plane that define each of the dif- 
ferent stellar populations are indicated with solid lines, as ded- 
scribed in the text. Absolute magnitudes are shown, converted 
from apparent m agnitude s using the distance modulii and extinc- 
tions provided in lDottei^ ct al. (2010). Labels for blue straggler, 
red giant branch, horizontal branch and main-sequence turn-off 
stars are indicated. Stars with large photometric errors have been 
omitted from this plot. 
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Figure 2. RA and Dec coordinates for all stars in the GC NGC 
6205. Circles corresponding to one, two, three and four core radii 
are shown. 



with high concentration parameters (|Harrislll99d ) for which 
King models are known to provide a poor fit, we multiphed 
the total stellar light by a mass-to-light ratio of 2 in or- 
der to obtain estimates for the total stellar mass contained 
within each circle. Calculating the total stellar mass con- 
tained within each circle from King models requires a num- 
ber of assumptions that we will discuss fully in Section |4] 



3 RESULTS 

In this section, we present our catalogue along with the 
results of our comparisons between the sizes of the differ- 
ent stellar populations and the total stellar mass contained 
within each circle and annulus. 

3.1 Catalogue 

The numbers of BS, RGB, HB and MSTO stars found within 
several different circles are shown for all clusters in Ta- 
ble [T] along with the total number of stars with magnitudes 
brighter than 0.5 mag below the MSTO. Number counts are 
only shown whenever the spatial coverage is complete within 
the indicated circle. 

3.2 Population Statistics 

How can we use our catalogue to learn which, if any, clus- 
ter properties affect the appearance of CMDs? One way to 
accomplish this is by plotting the size of a given stellar pop- 
ulation in a particular circle versus the total stellar mass 
contained within it. From this, lines of best fit can be found 
that provide equations relating the size of each stellar popu- 
lation to the total stellar mass contained within each circle. 
As described below, this is ideal for probing the effects of 
the cluster dynamics on the appearance of CMDs. 

The rate of two-body relaxation for a cluster can be 
appro ximated using the half-mass relaxation time l|Spitzeij 

ligsTh : 

Uh = 1.7 X W^[rh{pc)f^N'/''[m/MQr^^''years, (1) 

where is the half-mass radius, N is the total num- 
ber of stars within and m is the average stellar mass. 
The half-mass radii of MW GCs are remarkably similar 
independent of mass, and simulations have shown that 
changes by a factor of at most a few over th e course of 
a cluster's lifetime (|Murravl [l009l : lHenon|[l973h . The GCs 
that comprise our sample show a range of masses span- 
ning rou ghly 3 orders of magn itude, and have comparably 
old ages (jPe Angeli et ahlboOSl ). Therefore, Equation [1] sug- 
gests that the degree of dynamical evolution (due to two- 
body relaxation) experienced by a cluster is primarily de- 
termined by the total cluster mass for the GCs in our sam- 
ple. In particular, more massive clusters are less dynami- 
cally evolved, and vice versa. Consequently, if the size of a 
given stellar population is affected by two-body relaxation, 
the effects should be the most pronounced in the least mas- 
sive clusters in our sample. Additionally, the rate of (di- 
rect) steUaj_colhsi£nsJriCT cluster mass 
(e.g. iDavies. Piotto fc De Angeiill2004l '). This suggests that, 
if a given stellar population is affected by collisions, the ef- 
fects should be the most pronounced in the most massive 
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Figure 3. The logarithm of the number of stars belonging to 
each stellar population is shown for each circle as a function of 
the logarithm of the total stellar mass. From left to right and top 
to bottom, each frame corresponds to number counts contained 
within a circle having a radius of Vc, 2rc, 3rc and Arc- Blue cor- 
responds to blue stragglers, red to red giant branch stars, green 
to horizontal branch stars and black to main-sequence turn-off 
stars. Estimates for the total stellar mass contained within each 
circle were found using single-mass King models, as described in 
Section [231 

clusters in our sample. Therefore, by comparing the size of 
each stellar population to the total stellar mass contained 
within a given circle, the effects of the cluster dynamics can 
be quantified. This technique also ensures a normalized and 
consistent comparison since it accounts for cluster-to-cluster 
differences in the fractional area sampled by the ACS field 
of view. That is, we are consistently comparing the same 
structural area for each cluster. The validity and implica- 
tions of all of these assumptions will be discussed further in 
Section g] 

Plots showing the number of stars belonging to each 
stellar population as a function of the total stellar mass 
contained within each circle are shown in Figure [S] Uncer- 
tainties for the number of stars belonging to each stellar 
population were calculated using Poisson statistics. We also 
plot in Figure[4]the number of stars belonging to each stellar 
population as a function of the total stellar mass contained 
in each annulus outside the core. That is, we considered 
the populations for each annulus individually, as opposed 
to considering every star with a distance from the cluster 
centre smaller than the radius of the outer-most circles. Re- 
call that we have neglected clusters for which our theoretical 
King models provide a poor description of the true density 
distributions. This was the case for clusters in our sample 
having a high concentration pa ramete r , mos t of which are 
labelled as post-core collapse in iHarri j l| 19961 ). 

We performed a weighted least-squares fit for every re- 
lation in Figure [3] and Figure [l] Slopes and y- intercepts for 



Figure 4. The logarithm of the number of stars belonging to 
each stellar population is shown for each annulus as a function of 
the logarithm of the total stellar mass. The annulus and colour 
corresponding to each inset and stellar population, respectively, 
are the same as in Figure |3] 



these lines are shown in Table [2] and Table O respectively. 
Uncertainties for the slopes and y-intercepts were found us- 
ing a bootstrap methodology in which we generated 1,000 
fake data sets by randomly sampling (with replacement) 
number counts from the observations. We obtained lines of 
best fit for each fake data set, fit a Gaussian to the subse- 
quent distribution and extracted its standard deviation. 

As shown in Tabled the power-law index is sub-linear 
for BSs within the core at much better than the 3 — a con- 
fidence level, and it is consistent wi th the slope obtained 
in ou r earlier analysis presented in iKnigge. Leigh fc Sills! 
(|2009l ). The slopes for the BSs are also sub-linear at bet- 
ter than the 3 — a confidence level for all circles outside the 
core. This is also the case for all annuli outside the core, 
as shown in Table [S] Note, however, that the uncertainties 
for the BS slopes are very large for all annuli outside the 
core, whereas this is not always the case for corresponding 
circles outside the core. This is the result of the fact that 
the number of BSs drops off rapidly outside the core in sev- 
eral clusters so that the corresponding Poisson uncertainties, 
which are given by the square-root of the number of BSs, are 
significant. The rapid decline of BS numbers with increas- 
ing distance from the cluster centre in these clusters has also 
contributed to an increased degree of scatter in the relations 
for annuli outside the core relative to the corresponding re- 
lations for circles outside the core. 

The slopes are consistent with being linear for all other 
stellar populations in the core within their respective 3 — cr 
confidence intervals. Thi s agrees with the resul t s of o ur ear- 
lier analysis presented in lLeigh. Sills fc Kniggel (|2009l ) when 
we performed the comparison using the total core masses. 
The slopes are also consistent with being linear for all cir- 
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cles outside the core for both HB and MSTO stars. The 
power-law indices are sub-hnear at the 3 — a confidence 
level only for RGB stars, and this is only the case for circles 
outside the core. The power-law index is nearly unity for 
the core RGB population, yet the associated uncertainty is 
very large. Upon closer inspection, the distribution of power- 
law indices obtained from our bootstrap analysis for RGB 
stars in the core is strongly bi-modal, with comparably-sized 
peaks centred at ~ 0.82 and ~ 1.0. This bi-modality is 
most likely an artifact of our bootstrap analysis caused by a 
chance alignment of data points in the log Mcore-log N_rgb 
plane. Upon performing the comparison for only those stars 
found within particular annuli, our results suggest that the 
slopes are consistent with being sub-linear for all stellar pop- 
ulations at the 3 — a confidence level in only the annulus 
immediately outside the core (i.e. rc < r < 2rc). The slopes 
are consistent with being linear for RGB, HB and MSTO 
stars in all other annuli. 

We also tried performing the same comparisons using 
the total number of stars in each circle and annulus as a 
proxy for the total stellar mass. In this case, the slopes are 
sub-linear for BSs within all circles and annuli at the 3 — cr 
confidence level. Once again, the uncertainties are very large 
for all annuli outside the core, whereas this is not the case for 
corresponding circles outside the core. The slopes are con- 
sistent with being linear within the 1 — cr confidence inter- 
val for all other stellar populations in all circles and annuli. 
Our results are therefore inco nsistent with those presented in 
iLeigh. Sills fc Kniggel (120091 ') for the core RGB populations, 
in which we found that RGB numbers scale sub-linearly with 
the number of stars in the core at the 3 — cr confidence level. 
We will discuss the implications of these new results in Sec- 
tion g] 

3.3 Blue Stragglers and Single-Single Collisions 

As a check of o u r pr evious results reported in 
iKnigge. Leigh fc SiUsI (|2009l ). we also looked for a cor- 
relation between the observed number of BSs in the cluster 
core and the number predicted from single-single (1-1-1) 
collisions. The results of this comparison are shown in 
Figure [S] We define the predicted number of BSs formed 
from 1-1-1 collisions as A^'i+i = tbs/t\+i, where tbs is the 
average BS lifetime and ti+i is the average time between 
1-1-1 collisions in the cluste r core. We adopt th e sam e 
definition for ti+i as used in iKnigge. Leigh fc SiUsI l|2009l ). 
and assume tbs = 1-5 Gyrs as well as an average stellar 
mass and radius of 0.5 M© and 0.5 R©, respectively. We 
also adopt a constant mass-to-light ratio of M/L = 2 
for all clusters. Central lumi nosity d ensiti es an d velocity 
disper sions were taken from iHarrisI (Il996l ) and IWebbinkI 
l|l985h , respectively. 

Upon performing a weighted line of best fit for every 
cluster in our sample that overlaps with the catalogue of 
IWebbinkI l|l985l ), we find a power-law index of 0.15 ± 0.03 
(the uncertainty was found using the bootstrap methodology 
described in Section I3.2|l . For the subset of dense clusters 
having a central luminosity density satisfying log po > 4, we 
find a power-law index of 0.36 ± 0.14. As before, we find no 
significant correlation with collision rate, even for the subset 
of dense clusters. Although we do find a weak dependence 
of BS numbers on collision rate for the entire sample, this 
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Figure 5. The logarithm of the number of BSs predicted to have 
formed in the core from single-single collisions Ni+i versus the 
logarithm of the observed number of BSs in the core Nsg. Filled 
circles correspond to clusters having central luminosity densities 
satisfying log po < 4, whereas open triangles correspond to dense 
clusters for which log po > 4. The adopted definition for Ni+i 
has been provided in the text. 



is not unexpected since the collision rate and the core mass 
are themselves correlated, and our results suggest that there 
exists a strong correlation between BS numbers and the core 
masses. 



4 SUMMARY & DISCUSSION 

In this paper, we have presented a catalogue for BS, RGB, 
HB and MSTO stars in a sample of 35 GCs. Our cata- 
logue provides number counts for each stellar population 
within several different circles centred on the cluster cen- 
tre. The radii of the circles were taken to be integer mul- 
tiples of the core radius, and we have focussed on the in- 
ner four core radii since the field of view extends beyond 
this for only a handful of the clusters in our sample. Par- 
ticular consideration was given to our selection criteria for 
the different stellar populations in order to ensure that 
they were applied consistently from cluster-to-cluster. In 
particul ar, we have improved upon our previous selection 
criteria l|Leigh. Sills fc Kniggel |2007| ) by fitting theoretical 
isochrones to the cluster CMDs. This provides an unam- 
biguous definition for the location of the MSTO, which acts 
as the primary point of reference for the application of our 
selection criteria. As a result, our new catalogue is highly 
homogeneous. 

We have used our catalogue to quantify the dependence 
of the size of each stellar population on the total stellar mass 
enclosed within the same radius. As described in Section IX^ 
this provides a means of quantifying the effects, if any, had 
by the cluster dynamics in shaping the appearance of CMDs 
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above the MSTO. Below, we summarize the implications of 
our results for each of the different stellar populations. 

4.1 Blue Stragglers 

We have confirmed our previous result that the numbers 
of BSs i n the cores of GCs scale su b-linearly with the core 
masses ijKnigge. Leigh fc Sillsl [2OO9I ) . That is, we find pro- 
portionately larger BS populations in low-mass GCs. There 
exist several possibilities that could explain the origin of 
this sub-linear dependence. First, we previously suggested 
that this could be an artifact of an anticorrelation be- 
tween the binary fraction and the cluster (or core) mass 
(|Knigge. Leigh fc Silldlioogh . This assertion stems from the 
fact that, if BSs have a binary origin, we expect their 
numbers to scale with the core mass as Nbs ~ fbMcore, 
where /(, is the binary fraction in the core. As before, 
we find that Nbs ~ AfcorfT"'^- Our result could there- 
fore be e xplained if ^ Mc J^ re^'^' ^K Second, we sug- 
gested in iLeigh. Sills fc Kniggd l|2009l ) that the fact that 
the least massive GCs in our sample should be more dy- 
namically evolved than their more massive counterparts 
could be contributing to the observed sub-linearity for BSs. 
In particular, the low-mass clusters in our sample should 
have experienced a significant depletion of their very low- 
mass stars as a result of stellar evaporation induced by 



two-body relaxation (e.g. Soitzer 19871: 


Heggie & Hutll2003l: 


iDe Marchi. Paresce & Portegies ZwartI 


2OIOI). In turn, this 



could contribute to a higher fraction of merger products hav- 
ing masses that exceed that of the MSTO in low-mass GCs. 
As a result, more merger products would appear brighter 
and bluer than the MSTO in these clusters' CMDs, leading 
to more merger products being identified as BSs. Finally, 
mass segregation could also be contributing to the observed 
sub-linear dependence for BSs. Again, this is the result of 
the fact that the rate of two-body relaxation, and therefore 
dynamical friction, is in general the fastest in low-mass clus- 
ters. BSs are among the most massive stars in GCs, so they 
should rapidly migrate into the core via dynamical friction 
in clusters for which the half-mass relaxation time is shorter 
than the average BS lifetime. It follows that proportionately 
more BSs could have drifted into the core via dynamical 
friction in low-mass GCs. This would also contribute to the 
observed sub-linear dependence of BS numbers on the core 
masses. 

This last hypothesis can be tested by comparing our 
scaling relations for progressively larger circles outside the 
core. If mass segregation is indeed the cause of the observed 
sub-linear dependence of BS numbers on the core masses, 
then we might expect the power-law index to systematically 
increase as we consider progressively larger circles. That is, 
we could be including more BSs that have not yet migrated 
into the core via dynamical friction, particularly in the most 
massive clusters in our sample. However, our results suggest 
that the power-law index remains roughly constant for all 
circles. This is the case for both the comparison to the to- 
tal stellar masses as well as to the total number of stars 
contained in each circle. The fact that these relations are 
comparably sub-linear within all circles can be interpreted 
as evidence that mass segregation is not the dominant effect 
contributing to the observed sub-linear dependence of BS 
numbers on the total stellar masses (or the total number of 



stars). We note that for many of the clusters in our sample, 
the spatial coverage is comparable to or exceeds the half- 
mass radius. This is a sufficiently large fraction of the total 
cluster area for our comparison to be sensitive to the effects 
of mass segregation. 

On the other hand, several GCs are known to 
show evidence for a bi-modal BS radia l distribution (e.g. 
iMapelh et al] I2OO6I : iLanzoni et all 120071 '). That is, in these 
clusters the number of BSs is highest in the central cluster 
regions and decreeises with increasing distance from the clus- 
ter centre until a second rise in BS numbers occurs in the 
cluster outskirts. This seconary outer peak has been shown 
to occur at a distance from the cluster centre that exceeds 20 
core radii in several cases. Consequently, the spatial coverage 
provided by the ACS data likely does not extend sufficiently 
far in most clusters to detect any bi-modality in the BS 
radial distribution. Nonetheless, if applied to a statistically- 
significant sample for which the spatial coverage is complete 
out to the tidal radius, the technique we have presented in 
this paper could provide a powerful constraint for the origin 
of the bi-modal BS radial distribution observed in several 
MW GCs by addressing the role played by mass segrega- 
tion. 

We also tried to correlate the number of BSs observed 
in the cluster core wi th the number predicted fro m single- 
single collisions. As in lKnigge. Leigh fc Silld (|2009l ). we find 
that BS numbers depend strongly on core mass, but not 
on collision rate. This also proved to be the case for the 
subset of dense clusters satisfying log po > 4. Our previous 
interpretation that our results provide strong evidence for 
a binary, as opposed to a coUisional, origin for BSs in GCs 
therefore remains the same. 

4.2 Red Giant Branch Stars 

The technique used in this paper to compare the sizes 
of the different stellar pop ulations was first presented in 
iLeigh. Sills fc Kniggd (|2009l ). In that study, we introduced 
our method and applied it to a sample of 56 GCs taken from 
iPiotto et al.l l|2002il in order to study their RGB populations. 
Our results were consistent with a sub-linear dependence of 
RGB numbers on the core masses. In particular, we found 
evidence for a surplus of RGB stars relative to MSTO stars 
in the cores of low-mass GCs. We considered several possi- 
ble causes for this result, but concluded that our analysis 
should ideally be repeated with superior photometry in or- 
der to properly assess effects such as completeness. Given 
the high-quality of the ACS data, we are now in a position 
to reassess our previous result for RGB stars. 

Upon applying our technique to the ACS sample, we 
find that the numbers of RGB stars scale linearly with 
the core masses to within one standard deviation. This is 
also the case for our comparison to the total number of 
stars in the core. This suggests that we should reject our 
previous conclusions f or thi s stellar population reported in 
iLeigh. Sills fc Kniggd (Hooi). Specifically, if we take a strict 
3 — a limit as our criterion for whether or not the slopes 
are sub-linear at a statistically significant level, then the 
core RGB slope reported in this paper is consistent with be- 
ing linear wherea s this was not the cas e for t he core RGB 
slope reported in iLeigh. Sills fc Knig^ (|2009l ). However, if 
we take a more stringent criterion for statistical significance. 
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then there is no inconsistency between our new and old re- 
sults for RGB stars, and both slopes are consistent with 
being linear. 

We also tried comparing the RG B catalogue presented 
in this paper with the one presented in lLeigh. Sills fc Kniggd 
This showed that the old RGB numbers are slightly 
deficient relative to the new numbers at the high-mass end. 
Although this difference is not sufficiently large to com- 
pletely account for the difference in slopes (for the com- 
parison with the core mass) found between our new and old 
RGB catalogues, it works in the right direction and is likely 
a contributing factor. If our uncertainties are also factored 
in, then our new and old slopes agree to within one standard 
deviation (due mainly to the large uncertainty for the new 
slope). The source of the disagreement between our old and 
new RGB catalogues is unclear, and we cannot say whether 
or not incompleteness (in the old data set) is the culprit. The 
results of our artificial star tests have at least confirmed that 
incompleteness is not an issue for our new catalogue, how- 
ever it could certainl y have contributed to the lo wer RGB 
numbers reported in iLeigh. Sills Kniggd (j 20091 ). Indeed, 
the central cluster density tends to be higher in more mas- 
sive clusters, which should negatively affect completeness. 
It is not clear, however, why this would have affected RGB 
stars more than MSTO stars in the old data set. Given all 
of these considerations, we feel that our new results show 
that this issue needs to be looked at in more detail before 
any firm conclusions can be drawn. 

The evidence in favour of RGB numbers scaling linearly 
with the core masses is interesting. For one thing, it suggests 
that two-body relaxation does not significantly affect RGB 
population size relative to other stellar populations of com- 
parable mass in even the dense central regions of GCs. This 
is not surprising, since two-body relaxation is a long-range 
effect for which the stellar radius plays a negligible role. Sec- 
ond, it suggests that collisions do not significantly deplete 
RGB stars relative to other stellar populations despite their 
much larger radii. This is because the collision rate increases 
with increasing cluster mass, so we would expect RGB stars 
to appear preferentially depleted in massiv e clusters if they 
are s i gnificantly affected by collisions ( e.g. iBeers fc Davied 
12004 Ibavies. Piotto fc De Angelil 120041 '). Third, it suggests 
that the sub-linear relation found for BSs does not con- 
tribute to a sub-linear relation for RGB stars despite the 
fact that BSs should eventua lly evolve to occupy our RGB 
selection box, as discussed in lLeigh. Sills fc Kniggd (|2009l ). 
This is likely the result of the relatively small sizes of the BS 
populations in our sample when compared to the numbers 
of RGB stars since the rate at which evolved BSs ascend the 
RGB is thought to b e comparable to the RGB lifetim es of 
regular MSTO stars l|Sills. Karakas fc Lattanzidl2009l '). Al- 
ternatively, this result could, at least in part, be explained 
if a smaller fraction of BSs end up sufficiently bright and 
blue to be identified as BSs in the CMDs of massive GCs. 
In other words, it could be that a larger fraction of BSs 
are hidden along the MS in massive clusters, as discussed in 
Section [4. II In this case, the contributions to RGB popula- 
tions from evolved BSs could be comparable in all clusters, 
in which case a linear relationship between RGB numbers 
and the core masses would be expected. Finally, evolved BSs 
would be expected to have a negligible impact on RGB pop- 
ulation size if the average BS lifetime is considerably longer 



than the lifetimes of RGB stars. This effect is difficult to 
quantify, however, given tha t BS lifetimes are poorly con- 
strainod in the literature (e.g. lSandguist. Bolte fc HernquistI 
ll997. : .Sills et al..,2001i ) . 

4.3 Horizontal Branch Stars 

Our results suggest that HB numbers scale linearly with 
the core masses. This can be interpreted as evidence that 
two-body relaxation does not significantly affect the radial 
distributions of HB stars in GCs relative to the other stellar 
populations above the MSTO. One reason to perhaps expect 
that two-body relaxation should affect the spatial distribu- 
tions of HB stars stems from the fact that RGB stars are 
among the most massive stars in clusters, and they undergo 
significant mass loss upon evolving into HB stars. Conse- 
quently, the progenitors of HB stars should be heavily mass 
segregated. HB stars themselves, however, have relatively 
low-masses so that two-body relaxation and strong dynami- 
cal encounters should act to re-distribute them to wider or- 
bits within the cluster potential. The HB lifetime is roughly 
constant at 10® years (Ibcn ,19911) and it is comparable to 
or exceeds the core relaxatio n times for m ost of the low- 
mass clusters in our sample (|Harrisl [l996l ). Therefore, we 
might expect the HB populations in these clusters to ex- 
hibit more extended radial profiles relative to more massive 
clusters. This would contribute to a sub-linear relationship 
between the numbers of HB stars and the core masses. Our 
uncertainties are sufficiently large that this possibility can- 
not be entirely ruled out, however our results are consistent 
with a linear relationship between HB numbers and the core 
masses. 

4.4 Additional Considerations 

Recent observations have revealed the presence of mul- 
tiple stellar popula tions in a number of MW GCs (e.g. 
iPancino et al.|[2003l ). The majority of these cases have been 
reported in very massive clusters. Moreover, their existence 
is thought to be related to the chemical properties of GCs, in 
particular an observed anticorrelation between their sodium 
and oxygen abundances. In turn, these chemical signatures 
have bee n argued to be linked to the cluster metallicity, mass 
and age (Carr etta et al.l^2010^ . 

We identified clusters in our sample currently known 
to host multiple stellar populations, but none of these were 
clear outliers in our plots. Consequently, the effects had on 
our results by multiple stellar populations remains unclear. 
It is certainly possible that multiple stellar populations have 
contributed to the uncertainties for the weighted lines of best 
fit performed for the relations in Figure [3] and Figure |31 It 
is difficult to quantify the possible severity of this effect, 
however, given the limited evidence linking multiple stellar 
populations to cluster properties. 

Although the uncertainties are sufficiently large that the 
slopes are consistent with being linear at the 3 — a confidence 
level for all stellar populations when performing the compar- 
isons with the total stellar mass, the reported slopes are typi- 
cally less than unity within the 1— cr, and often even the 2— a, 
confidence interval. This does not appear to result from the 
fact that we have obtained our estimates for the total stel- 
lar masses by numerically integrating 3-dimensional density 
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distributions and are comparing to number counts, which 
are projected quantities. To address this, we also tried ob- 
taining the total stellar masses by numerically integrating 2- 
dimensional surface brightness profiles so that we are consis- 
tently comparing only projected quantities. Despite this, our 
results remain unchanged and the new slopes agree with the 
old ones to within one standard deviation for all stellar pop- 
ulations. Another possibility to account for this trend that is 
perhaps worth considering is a systematic dependence of the 
mass-to- light ratios of clusters on their total mass. There are 
two ways this could have affected our analysis. First, stellar 
remnants have been shown to affect the dynamical evolu- 
tion of clusters, an d therefor e the sizes of their cores (e.g. 
Lee. Fahlman fc Richer 199ll : iTrenti. Vesperini fc Pasguatd 
It follows that, if the number of stellar remnants de- 
pends on the cluster mass, then this could contribute to an 
additional underlying dependence of the core radius on the 
cluster mass. This could perhaps arise as a result of the fact 
that the ratio of the rate of stellar evolution to the rate of 
dynamical evolution is larger in more massive clusters, since 
in general the rate of two-body relaxation decreases with 
increasing cluster mass whereas the rate of stellar evolu- 
tion is independent of the cluster mass. Coupled with their 
deeper gravitational potential wells, this could contribute 
to more massive clusters retaining more stellar remnants. 
Second, variations in the mass-to-light ratios of clusters can 
also occur as a result of changes in the average stellar mas s 
(not including stellar remnants) (|Kruiissen fc MieskdIioogI ). 
That is, we can approximate the total stellar mass contained 
in the core as: 



Mcor 



4 M 3 



mNcor 



(2) 



where M/L is the mass-to-light ratio, po is the central lu- 
minosity density, m is the average stellar mass and Ncore is 
the total number of stars in the core. Based on our results, 
Mcore oc N^ore: whcrc wc have used the total number of stars 
in the core with magnitudes brighter than 0.5 mag below the 
MSTO as a proxy for Ncore- This could suggest that m oc 
^core OC M~Jl,} . In othcr words, the average stellar mass in 
the core decreases weakly with increasing core mass. This 
could in part be due to the fact that more massive clusters 
should be less dynamically evolved than their less massive 
counterparts, and should therefore be less depleted of their 
low-mass stars due t o stellar eva poration induced by two- 
body relaxation (e. g. Ambartsumian.l938i: Spitzer fc HarmI 



19581 : iHenonI Il960l : iDe Marchi. Paresce fc Portegies ZwartI 



20101 ). Similarly, mass segregation should also tend to op- 



grate preferentially massive stars into the core ( 


SDitzerlll96g|: 


Spitzer fc Hard 1971 : Farouki fc Salpeter 1982 


: IShara et al.l 


1995l:lKing. Sosin fc Coolll 19951: iMevlan fc Heegielll997l). Al- 



ternatively, differences in the stellar mass function in the 
core could result from variations in the degree of primor- 
dial mass segregation, or even variations in the initial stellar 
mass function. 

We have assumed throughout our analysis that the core 
mass is a suitable proxy for the total cluster mass. We have 
checked that these two quantities are indeed correlated, how- 
ever this does not tell the whole story since we are also using 
the total cluster mass as a proxy for the degree of dynamical 
evolution. The central concentration parameter, defined as 
the logarithm of the ratio of the tidal to core radii, describes 



the degree to which a cluster is centrally concentrated. Pre- 
vious studies have shown that there exists a weak correlation 
between t he concentrati on parameter and the total clu ster 
mass (e.g. lDiorgovski fc Mevlan.l994. : ,McLaughlinll2000l '). In 
order to better use our technique to reliably probe the effects 
of the cluster dynamics on the sizes and radial distributions 
of the different stellar populations, the concentration pa- 
rameter should ideally be accounted for when applying our 
normalization technique in future studies. It is not yet clear 
how the concentration parameter can be properly absorbed 
into the normalization, however its effect on our analysis 
should be small given the weak dependence on cluster mass. 

The assumption that the degree of dynamical evo- 
lution experienced by a given cluster depends only on 
its mass is also incorrect. Two-body relaxation has been 
shown to dominate cluster evolution for a significa nt frac- 
tion o f the lives of old MW GCs (e.g. iGielos. Hcggie fc Zhaol 
l201ll ). however other effects can also play a significant 
role. For example, stellar evolution is known to affect 
the dynamical evolution of star clusters, although its 
primary role is played dur i ng their early ev olutionary 
phases (e.g. lApplegatd Il986l : IChernoff fc Weinb erg 1990|; 
iFukushige fc Heggid Il995l ). Tidal effects from the Galaxy 
have also been shown to play an important role in decid- 
ing the dynamical fates of clusters by incre asing the rate 
of mass loss across the tidal boundary (e.g. iHeggie fc HutI 
Consequently, clusters with small perigalacticon dis- 
tances should appear more dynamically evolved than their 
total mass alone would suggest. This effect can be signif- 
icant, and has likely contributed to increasing the uncer- 
tainties found for the comparisons to the total stellar mass. 
Therefore, tidal effects from the Galaxy should also ideally 
be absorbed into our normalization technique in future stud- 
ies. This can be done by using the perigalacticon distances 
of clusters as a rough proxy for the degree to which tides 
from the Gala xy should have affected their internal dynam- 
ical evolution jCieles. Heggie fc Zhad[201ll ). 

Interestingly, tides could also help to explain why the 
uncertainties for the comparisons to the total number of 
stars in each circle are considerably smaller than for the 
comparisons to the total stellar mass. We have used the 
number of stars with magnitudes brighter than 0.5 mag be- 
low the MSTO as a proxy for the total number of stars. 
Consequently, we are comparing stars within a very narrow 
mass range, so that all populations of interest should have 
been comparably affected by two-body relaxation (except, 
perhaps, for HB stars) independent of tidal effects from the 
Galaxy. In other words, tides should affect all stars above 
the MSTO more or less equally, and this is consistent with 
our results. It is also worth mentioning here that our King 
models consider only a single stellar mass. This assumption 
is not strictly true and could also be contributing to increas- 
ing the uncertainties found for the comparisons to the total 
stellar masses. 

An additional concern is that we do not know if the 
clusters in our sample are currently in a phase of core con- 
traction or expansion. This has a direct bearing on the re- 
cent history of the stellar density in the core, and therefore 
the degree to which stars in the core should have been af- 
fected by close dynamical interactions. These effects are in- 
dependent of two-body relaxation and occur on a time-scale 
that is typically much shorter than the half-mass relaxation 
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time (|Heggie fc HutI 120031 '). The effects could be significant 
in clusters that were recently in a phase of core-collapse 
but have since rebounded back out of this highly concen- 
trated state. This could occur, for example, as a result of 
binary formation induced by 3-body interactions combined 
with their subsequent hardening via addit ional encounters 
l|Hut fc Bahcalllll983l : iHeggie fc Hutll2003l ). In general, bi- 
naries play an important role in the dynamical evolution of 
clusters, and could have affected our results in a number of 
ways. This is a difficult issue to address even qualitatively 
given how little is currently known about the binary popula- 
tions in globular clusters. Theoretical models suggest, how- 
ever, that the time-scale for core contraction is often longer 
than a Hubble time, and that this evolutionary phase will 
only come to an end once the central density becomes suf- 
ficiently high for hardening encounters involving binaries to 
halt the process (Frcgcau, Ivanova fc Rasio 2009). It follows 
that the cores of most MW GCs are expected to currently 
be in a phase of core contraction. This process is ultimately 
driven by two-body relaxation, so that our assumption that 
the total cluster mass provides a suitable proxy for the de- 
gree of dynamical evolution is still valid. 

In summary, our results suggest that effects related to 
the cluster dynamics do not significantly affect the relative 
sizes of the different stellar populations above the MSTO. 
This is the case for at least RGB, HE and MSTO stars. 
BSs, on the other hand, show evidence for a sub-linear de- 
pendence of population size on the total stellar mass con- 
tained within the same radius. Whether or not the clus- 
ter dynamics is responsible for this sub-linearity is still not 
clear. Notwithstanding, our results have provided evidence 
that mass segregation is not the dominant cause for this re- 
sult, although it will be necessary to redo the comparison 
performed in this study with a larger spatial coverage in or- 
der to fully address this question. Further insight into the 
origin of the sub-linearity found for BSs will be provided by 
reliable binary fractions for the clusters in our sample, which 
are forthcoming (Sarajedini 2010, private communication). 
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Table 1. Stellar Population Catalogue 
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Table 2. Lines of Best Fit for log(Mj,j^c(e/10^) Versus log(Npop) 



Circle BS 



RGB 



HB 



MSTO 



<re log(Nss) = (0.39 ± 

0.05)log(M<rJ103) + 
(1.22 ± 0.05) 

< 2rc log(Nss) = (0.36 ± 

0.05)log(M<2rJ10^) 

+ (1.26 ± 0.08) 

< arc log(NBs) = (0.47 ± 

0.04)log(M<3^J103) 
+ (1.02 ± 0.08) 

< 4rc log(NBs) = (0.45 ± 

0.05)log(M<4rJ103) 
+ (1.01 ± 0.12) 



log(NHGB) = (0.95 ± 

0.11)log(M<rJ103) + 
(1.36 ± 0.11) 

log(NflGB) = (0.87 ± 

0.08)log(M<2r,/10=^) 

+ (1.27 ± 0.12) 
log(NHGB) = (0.80 ± 
0.06)log(M<3^J103) 
+ (1.26 ± 0.12) 
log(NHGs) = (0.75 ± 

0.07)log(M<4r,/10=') 

+ (1.28 ± 0.15) 



log(Nj^s) = (0.95 ± 
0.06)log(M<rJ103) + 
(0.92 ± 0.06) 
log(NHB) = (0.85 ± 
0.06)log(M<2rJ10:^) 
+ (0.84 ± 0.09) 
log(NffB) = (0.79 ± 
0.08)log(M<3^J103) 
+ (0.83 ± 0.14) 
log(NjfB) = (0.75 ± 
0.09)log(M<4rJ103) 
+ (0.83 ± 0.19) 



log(NMSTo) = (0.90 ± 
0.07)log(M<^J103) + 
(2.03 ± 0.07) 
log(NMSTo) = (0.82 ± 

0.06)log(M<2r,/103) 

+ (1.98 ± 0.09) 
log(NMSTo) = (0.82± 
0.11)log(M<3^J103) 
+ (1.86 ± 0.20) 
log(NMSTo) = (0.78± 
0.12)log(M<4^J103) 
+ (1.86 ± 0.25) 



Table 3. Lines of Best Fit for log{Mannulus/^0^) Versus 
log(Npop) 



Annulus 



BS 



RGB 



HB 



MSTO 



< r < 2rc 



2rc < r < 3rc 



3rc < r < 4rc 



log(NBs) = (0.27 ± 
0.08)log(M,^<,<2,J103) 
+ (1.04 ± 0.13) 
log(NBs) = (0.39 ± 

0.09)log(M2r.e <r.<3r., /lO^) 

+ (0.52 ± 0.16) 
log(NBs) = (0.15 ± 
0.21)log(M3r 

c <r<4rc 

+ (0.77 ± 0.35) 



log(N«GB) = (0.80 ± 
0.06)log(M,^<,<2,,/103) 
+ (1.22 ± 0.08) 
log(NflGs) = (0.79 ± 

0.11)log(M2r-,<r.<3r-,/103) 

+ (0.97 ± 0.17) 

log(NflGs) = (0.63 ± 
0.11)log(M3, 

c<r<4:rc 

/103) 

+ (1.04 ± 0.18) 



log(N^jB) = (0.77 ± 
0.06)log(M,^<,<2.,/10'^) 
+ (0.79 ± 0.09) 
log(NHB) = (0.68 ± 

0.10)log(M2r.,<r.<3re/10^) 

+ (0.67 ± 0.16) 
log(N^fB) = (0.69 ± 
0.18)log(M3, /103) 
+ (0.45 ± 0.31) 



log(NMSTo) = (0.76 ± 

0.05)log(M,,<,<2r,/10») 

+ (1.91 ± 0.08) 
log(NMSTo) = (0.78 ± 

0.12)log(M2r,<r-<3re/103) 

+ (1.64 ± 0.20) 
log(NMSTo) = (0.69 ± 

0.16)log(M3r,<r<4rJ103) 

+ (1.60 ± 0.27) 
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APPENDIX A: STELLAR POPULATION 
SELECTION CRITERIA 

In this section, we present our selection criteria for BS, RGB, 
HB and MSTO stars. Our met h od is similar to that de- 
scribed in iLeigh. Sills fc Kniggel ll2007t ). and we have used 
this as a basis for the selection criteria presented in this 
paper. First, we define a location for the MSTO in the 
(F606W-F8I4W)-F814W plane using our isochrone fits. The 
MSTO is chosen to be the bluest point along the MS of each 
isochrone, which we denote by {{y-l)MSTO,lMSTo)- In or- 
der to distinguish BSs from MSTO stars, we impose the 
conditions: 



F8UW s; mi {F606W - F8UW) + 611, (AI) 

where the slope of this line is mi — —9 and its y-intercept 
is given by: 

611 = (Imsto - 0.10) - mi((V - I)msto ~ 0.10) (A2) 

Similarly, we distinguish BSs from HB stars by defining 
the following additional boundaries: 



F8UW ^ 

F8UW ^ 

F8UW ^ 

F8UW ^ 

{F606W - F814:W) ^ 

F8UW ^ 



mi(F606W - F8UW) + 6i<A3) 
m2(F606W - F8UW) + 62i(A4) 
m2(F606W - F8UW) + 62<A5) 
m/fs (F606W - ^^8141^) + ^M) 
{V - I)hb (A7) 
Imsto, (A8) 
I)„B ^ {V - 



where m2 = 6, mns ~ —1.5 and {V 
I)msto — 0.4. We also define: 

612 = (/m5to -0.55) -mi ((1/-/)msto -0.55) (A9) 

621 = (Imsto- 0.80) -m2((l/-/)MSTO + 0.10) (AlO) 

622 = {Imsto + 0.30) - m2{(V - I)msto - 0.20),(A11) 

and bHB ~ Ihb + 1.2, where Ihb roughly corresponds to 
the mid-point of points that populate the HB and is chosen 
by eye for each cluster so that our selection criteria best fits 
the HB in all of the CMDs in our sample. 

We apply a similar set of conditions to the RGB in order 
to select stars belonging to this stellar population. These 
boundary conditions are: 



F814W^ ^ mHB(Fm&W ~ F8im) + bHB (A12) 

F814W ^ mHGs(F606W-J'814W)-|-&3i (A13) 

F814W^ mRGB{Fmm ~F8l'm) + hz2 (A14) 

F814iy ^ Ibgb, (A15) 

where itirgb = —23, Irgb is defined as the F814W magni- 
tude corresponding to a core helium mass of 0.08 M© and: 

631 = (Imsto -0.60) -mB.GB{{V-I)MSTO+0.i)S^16) 

632 = (Imsto -0.60) -mB.GB{{V~I)MSTO+0.tS:)17) 

Core helium-burning stars, which we refer to as HB 
stars, are selected if they satisfy one of the following sets 
of criteria: 



F8UW > 
F8UW s; 
(F606W - F814W) < 



mHBiF606W 
Imsto + 2.5 
{V — I)msto - 



- F814W) MMh) 
(A22) 
0.4, (A23) 



F814W < mi(F606Vy-F814W)-f &{A24) 
i^814W > mHB{F606W - F8im)M.Mb) 
(F606W - 7^814^) ^ {V - I)msto - 0.4 (A26) 

We define (V — /)hs on a cluster- by-cluster basis in order 
to ensure that we do not over- or under-count the number 
of HB stars. This is because the precise value of (F606W- 
F814W) at which the HB becomes the RGB varies from 
cluster-to-cluster. In addition, the precise location of the 
transition in the cluster CMD between HB and EHB stars 
remains poorly understood. To avoid this ambiguity, we con- 
sider HB and EHB stars together throughout our analysis, 
and collectively refer to all core helium-burning stars as HB 
stars throughout this paper. 

Finally, MSTO stars are selected according to the fol- 
lowing criteria: 

F814M/ > Irgb (A27) 
F814M/ > mx(F606W -F8l4W) + hix (A28) 
F8im ^ (V -I) msto (A29) 
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